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Planck scale physics may influence the evolution of cosmological fluctuations in the early stages 
of cosmological evolution. Because of the quasi-exponential redshifting, which occurs during an 
inflationary period, the physical wavelengths of comoving scales that correspond to the present 
large-scale structure of the Universe were smaller than the Planck length in the early stages of the 
inflationary period. This trans-Planckian effect was studied before using toy models. The Horava- 
Lifshitz (HL) theory offers the chance to study this problem in a candidate UV complete theory of 
gravity. In this paper we study the evolution of cosmological perturbations according to HL gravity 
assuming that matter gives rise to an inflationary background. As is usually done in inflationary 
cosmology, we assume that the fluctuations originate in their minimum energy state. In the trans- 
Planckian region the fluctuations obey a non-linear dispersion relation of Corley-Jacobson type. In 
the "healthy extension" of HL gravity there is an extra degree of freedom which plays an important 
role in the UV region but decouples in the IR, and which influences the cosmological perturbations. 
We find that in spite of these important changes compared to the usual description, the overall 
scale-invariance of the power spectrum of cosmological perturbations is recovered. However, we 
obtain oscillations in the spectrum as a function of wavenumber with a relative amplitude of order 
unity and with an effective frequency which scales nonlinearly with wavenumber. Taking the usual 
inflationary parameters we find that the frequency of the oscillations is so large as to render the 
effect difficult to observe. 

PACS numbers: 



I. INTRODUCTION 

The inflationary scenario [l[ has become the current paradigm to describe the early evolution of the universe. 
The success of this theory comes from a number of predictions which have been confirmed by the observations, the 
main one being the causal mechanism for the generation of the primordial cosmological perturbation which inflation 
provides @]. Inflation is generally modelled by General Relativity (GR) coupled to scalar field matter. The scalar 
field (called inflaton) is responsible for the accelerated expansion of space. The success of inflation as a theory for 
the origin of structure comes from the fact that the wavelengths of fixed comoving modes which are of cosmological 
interest today, namely wavelengths corresponding to the present large-scale structure and to the cosmic microwave 
background anisotropies, were sub-Hubble at the beginning of the period of inflation. They cross the Hubble horizon 
about 50 - 60 Hubble times, or e-folds, before the end of inflation. Thus, inflation must last at least that long to solve 
the problems of the standard cosmological model (SCM) [H, H|. 

Most inflationary models have a period of accelerated expansion which lasts much more than the 50 - 60 e-folds 
required solve the SCM problems. If the period of inflation last for more than about 70 e-folds (this number depends 
very slightly on the energy scale at which inflation takes place - to get the above number we are assuming that it is 
the Grand Unification scale) , then all of the wavelengths important today were smaller that the Planck length at the 
beginning of inflation. In this regime we cannot trust the calculations since General Relativity will no longer yield 
a good description of the physics. This is the known trans-Planckian problem for cosmological fluctuations Q. A 
correct computation of both the generation and evolution of fluctuations in the trans-Planckian regime must be done 
in the correct UV completion of the theory. 

This problem is named after the analogous trans-Planckian problem for Hawking radiation in black-hole physics 
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Q . In the cosmological context this problem was discussed in the literature in many different toy model modifications 
of General Relativity coupled to a canonical scalar field. These models describe the evolution on sub-Planckian 
wavelengths by imposing ad-hoc modifications of the usual physics, e.g. modified dispersion relations 0, H| or non- 
commutativity |j[ (see also [Hj). It is possible to construct toy models in which the standard predictions of inflation 
are changed, in others they are maintained Q either fully or up to small corrections. For example, imposing the usual 
initial conditions on a time-like "new physics hypersurface" given by the physical wavelength being some fixed 
UV cutoff scale smaller or equal to the Planck scale yields a scale-invariant power spectrum of curvature fluctuations 
with small amplitude superimposed oscillations 42j If, on the other hand, we assume that a short period of inflation 
is preceded by a short non-singular bounce and (before that) the time reverse of the background evolution, then a 
deep red spectrum with index n s = —3 results [14| . However, to go beyond toy model studies, we must be able to 
perform the analysis of cosmological fluctuations in a full quantum gravity theory. 

Hofava-Lifshitz (HL) gravity [15| is a candidate theory of quantum gravity which provides a well motivated frame- 
work for describing the evolution of the universe on trans-Planckian scales. It is an theory of gravity in 3+1 dimensions 
which uses the same metric degrees of freedom as General Relativity but which is power- counting renormalizable (with 
respect to the scaling symmetry to be introduced below) . This is achieved by picking a preferred time direction, thus 
abandoning the full space-time diffcomorphism invariance (reducing the symmetry to simply spatial diffcomorphisms), 
and by introducing an anisotropic scaling of Lifshitz type. The loss of space-time isotropy has as consequence the 
absence of Lorentz symmetry, a symmetry which appears as an emergent symmetry in the infrared (IR), where GR 
is recovered. 

As a consequence of the loss of the space-dependent time reparametrization symmetry of General relativity, one of 
the four gauge symmetries for cosmological perturbations (namely the space dependent time rescaling symmetry) is 
lost, and since the theory presents the same number of basic degrees of freedom as General Relativity, there is an 
extra scalar perturbation mode (scalar in terms of its transformation under spatial rotations) (l6j . 

Since it was proposed by Horava in its most simple version fl5j , many other versions of the theory were proposed 
in the literature |43|. We have mainly two versions: the "projectable" , where the lapse function depends only on time 
to mimic the reduced symmetry of the problem; or the " non-projectable" version, where the lapse function depends 
on space and time as in GR. In the " projectable" version the extra degree of freedom of the theory is dynamical (l8j 
and ghost-like or tachyonic. This version is also plagued by "strong coupling" problems fl9l. [20|. which may spoil the 
validity of this setup. 



In the original form of the "non-projectable" version, the extra degree of freedom is non-dynamical 21 1 at linear 
order about a spatially homogeneous and isotropic background. "Strong coupling" problems may arise 2(f [HI- The 
most studied type of a "non-projectable" theory (one which contains all of the terms which are allowed by power- 
counting renormalizability and symmetry) is the "healthy extension" [23| that appears to fix these "strong couplings" 
problems and pathologies associated with the extra degree of freedom that appear in the other versions of the theory 



[44( . Fluctuations in the healthy extension of HL gravity were first considered in [25[ . As it was shown in [26[ , the 
extra degree of freedom now becomes dynamical. With appropriate choices of the parameters in the HL Lagrangian, 
the extra degree of freedom is neither ghost-like nor tachyonic. It is important only in the UV and decouples in the 
IR. This is interesting from the point of view of cosmology since we know, from observations, that we can only have 
one propagating scalar metric degree of freedom in the IR. On the other hand, the fact that in the UV there are 
two propagating degrees of freedom makes the theory a promising one to study in the context of the trans-Planckian 
problem for fluctuations (4f|. 

In the presence of spatial curvature, HL gravity can lead to bouncing cosmological solutions [29j . Here, however, we 
will assume the existence of scalar field matter which leads to an inflationary phase during the evolution of the early 
universe. We will use the "healthy extension" of the HL gravity to study the trans-Planckian problem for fluctuations. 
We are interested in checking, based on a well motivated description of the UV and IR theory, if the predictions of 
inflation for cosmological perturbations are altered. We will have to take into account two effects: the change in the 
physics in the UV region that will be manifest in a modified dispersion relation of Corley-Jacobson type [3CJ (as the 
ones which were studied in Q); and the presence of an extra degree of freedom that will be dynamical in the UV (for 
wavelengths smaller that Hubble radius in our case). This extra degree of freedom can be interpreted as an "entropy" 
perturbation, in analogy with what occurs in multi-field inflation models [3lj. but in this case the extra field is a 
scalar perturbation of the metric. The entropy mode is coupled to the "adiabatic perturbation" and will influence its 
dynamics while the wavelength is sub-Hubblc. Thus, we need to analyse the coupled system of differential equations 
for both scalar perturbations modes in the different regimes of the theory in order to evaluate the final power spectrum 
of the " adiabatic" mode. 

The article is organized as follows. In Section|Hl we give a brief description of the "healthy extension" of HL, setting 
up the theory we are using and calculating linear perturbations about a homogeneous and isotropic background. We 
show how the reduced symmetry group gauges only one of the degrees of freedom, leaving one extra degree of freedom, 
by calculating the second order action. In Section |HT] we study the cosmological perturbations. First we find the 
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canonical variables of the problem and show that one of these corresponds to the usual Mukhanov-Sasaki variable and 
the other decouples in the IR. We take the UV limit of the theory to find the contributions of the trans-Planckian 
physics. The initial state is calculated by minimizing the vacuum energy and the solutions to the inhomogeneous 
equations of motion are found in each region of the problem. Finally, the power spectrum is calculated. We find 
that the overall scale-invariance of the spectrum is maintained, but that there are superimposed oscillations with an 
effective frequency depending on the wavenumber (very different from the oscillations which are seen in the approach 
of [ll| to the trans-Planckian problem). The amplitude of these oscillations is not suppressed (again in contrast 
to what is seen in However, the frequency of oscillation is very large making the effect extremely difficult to 

measure. 



II. HEALTH EXTENSION OF HORAVA-LIFSHITZ 



A. Background 

In HL theory it is natural to use the ADM (Arnowitt, Desner and Misner) formalism that separates space-time into 
time and spatial foliations: 

ds 2 = -N 2 dt 2 + gij (dx l - N l dt) (dx j - N j dt) , (1) 

where <?y (t, x) is the metric of the spatial section, N is the lapse function and N l (t, x) is the shift vector. As 
mentioned, the lapse function can be restricted to depend only on time, yielding the "projectable" theory, or else 
it is taken to depend on both space and time, yielding the "non-projectable" version. The "healthy extension" is a 
non-projcctablc version of the theory with all of the terms allowed by the residual symmetries and by power-counting 
renormalizability included. Thus N(t,x) and we describe the same lapse function as in GR. In this case, we must 
include in the action terms proportional to the quantity: 

diN 

* = — ■ (2) 

Because of the anisotropic scaling, the classical scaling dimensions of the coordinates of the theory (when the scaling 
dimension is z = 3 as it is in four space-time-dimensional HL gravity) are, in units of spatial momenta: [t] s = —3, 
1. With this we have that [V = did 1 } = 2. The ADM fields have scaling dimension: [gij] s = [N] s = and 



ijs — vviun unio wu nave luiat [v — u t u j — ^. iiib r-vj^ivi. iicilio nave oecimig, umiciioiuii. \y t j 

[Ni}=2. 

The action of HL gravity is constructed by including all terms with scaling dimension smaller or equal to six, in order 
to obtain a theory which is power-counting rcnormalizablc, and which arc invariant under the spatial diffcomorphisms, 
i.e. under elements of the symmetry group Diffjr. The "healthy extension" action can be written as, following (26j 
(see also HHH): 

S = X 2 J dtd 3 x^gN (£ km ~ C v ~ C E + X" 2 £m) , (3) 

where x 2 = l/lGirG. The kinetic term is the one that contains time derivatives of the metric in a covariant way with 
respect to the new foliating preserving symmetry, Diffj^: 

C km = K l3 K ij - \K 2 , (4) 

where Kij = (1/2N) (g^ — ViNj — VjJV,) is the extrinsic curvature of the constant time hypcrsurfaces. The potential 
term is constructed by including all terms that preserve Diffj? until the sixth spatial derivative: 

C v = 2A - R+ -j (g 2 R 2 + g 3 RijR ij ) + — [g 4 R 3 + g 5 RRijR ij + g&R)R{R k i 

+ g 7 RV 2 R + g s (V. t R jk ) (V l W k )] , (5) 

where A is the cosmological constant and has dimension of momentum squared and the coefficients are dimensionless. 
The new components proportional to cij that appear in the "healthy extension" are: 

C E = -va.^ + \ (r/ 2 a 4 Aa ?; + ?] 3 RV l a t + ...) + \ [^A V + rjsARV^ + r/ 6 i? 2 V 4 a l + ...], (6) 
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where the coefficients are dimensionless. As we are interested in obtaining an inflationary background cosmology, we 
include matter in the form of a scalar field, with matter Lagrangian [34j 

where the potential terms are given by: 

V (<f>, d^, 9ij ) = V + VxPq + V 2 Pf + V 3 P? + V^P 2 + V b P a P2 + VePiP 2 , (8) 

with Po = (V0) 2 , Pi = A l <ft and A = g^'Vi'Vj. The scalar field has scaling dimension [<fr] s = (d — z)/2 = 0. In the 
way it is written above, the coefficients are not dimensionless as the previous coefficients are. If we want to write 
them as dimensionless coefficients we have to make the following substitution: the order 2 coefficient V\ is already 
dimensionless; the order four coefficients are changed to V 2 — > V 2 /x 2 an d ^4 ~~ > ^i/x 2 ; with V 2 and V4 dimensionless; 
and the dimension six coefficients are changed to V3 — > V3/X 4 , V5 —> Vs/x 4 an d V§ —> Vq/x 4 , with V3, V5 and Vq 
dimensionless. 

The non-projcctable version of HL can contain more than 60 terms. Here we only considered the ones that are 
important for the scalar part of the perturbations, the perturbations we interested in. In the IR the GR action is 
recovered when A = 1. 

We can now obtain the Hamiltonian and momentum constraints. Varying the action with respect to N" 1 (t, a;- 7 ), we 
have the momentum constraint: 



K 2 



V <7 r« = — P , (9) 



where n ij = ^§P^ = -K lj + XKg ij and J { = -Nj$ = i - Nk ^k<f>) V 4 0. Note that © is a local constraint 



The Hamiltonian constraint can be obtained by varying the action with respect to the lapse function, N(t,x l ): 

S£ 1 

C cm + C pot +C E + N— y = —j J*, (10) 



SN 4 X 

with J* = 2 [N^- + C,fA and the new terms from "healthy extension": 



JS-— = 2r}\/ia ^-AVid H 5- A it r A Via H — jA KH s-AP + ... . (11) 



This constraint is a local constraint, unlike what happens in the "projectable" version of the theory. 

One might think that the presence of two local constraints should lead to a reduction of the number of degrees of 
freedom of the theory to have the same number as in General Relativity. However, this is not the case. The local 
Hamiltonian constraint is not able to decouple the extra degree of freedom in all the regimes of the theory, only in 
the IR. We can see this by studying cosmological perturbation theory. 

B. Perturbations 

Perturbing about an isotropic and homogeneous universe described by the spatially flat FRW metric, 

ds 2 = -dt 2 + ,r l \,,l.i-,l.v' , (12) 

the background values and scalar perturbation of the ADM variables are given by [||: 

N + SN (t, x k ) = 1 + v (t, x k ) , (13) 
N % + SNi (t, x k ) = d t B (i, x k ) , (14) 

g l3 + 5g l3 (t, x k ) = a 2 S lJ + a 2 (t) [-2^ (t, x k ) + 2P (t, x k ) .J . (15) 
We must also consider the perturbation in the scalar field, given by: 

<j>(t, x k ) = <j> (t)+5<j>(t, x k ) . (16) 
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The invariance under Dif 7>, allow us to set E = 0. Because of the reduced symmetry we can no longer constrain 
another degree of freedom of the perturbations and use the standard gauges used in GR 0]. Thus we are left with 
more degrees of freedom than we would have in GR: v, B, ip and 8<j). We can reduce the number of them by solving 
the constraints of the theory. 

Expanding the momentum and Hamiltonian constraints to first order, we can use them to solve only for two degrees 
of freedom (in a spatially flat background), namely v and B. In the momentum space, we can rewrite the equation 
using the physical momentum fc = k/a (as done in [26j ) : 



d(k) B k (t) = - (3A - 1) 
-|(3A-1) 



+ 2/i (k) 



4>k (f) - (3A - 1) 



X 2 k 2 



l k (t) - 4 (3A - 1) f 2 (fc) (t) 



Voj (<h) + Vi (0o) ~k 4 H + 3 (3A - 1) o ff 2 



§r " Mi (£) k 2 



5<t>k (t) 
X 2 k 2 



;i7) 



d (fc) v h (t) = (A - 1) ^5fa (t) - 4 (3A - 1) H^ k (t) + 4 (A - 1) / a (k) fcV* (t) + 



{(3A - 1) faH + (A - 1) [V 0t4 , (M + Vi (0 O ) fc 4 ] } 



5<t> k (t) 



(18) 



where we introduced the following functions to simplify these equations: 



k 2 fc 4 - k 2 fc 4 

h (k) = -v + m-2 + m-i , h (A) = -i + m-2 + m-r , 

A A A A 



cf (fc) = 4(3A- 1)H 
= 4(3A 



(A 



l)tf 2 



1) 
A 



^+2/i (A) P 

1 00^2 



2(3A^1) X 



A - 1 



2(3A- 1) 



fc 2 A: 4 

-f] + ??2^ + ?74^- 

X X 



H 2 



(19) 
(20) 
(21) 



Since we have the function d (fc) in both solutions, B k (t) and v k (t) are regular in the limit A — > 1 as long as H ^ 0. 
The second expression in (f2~Tj) is valid only when A ^ 1/3 and H ^ 0. From this, we can see that the value of 
the physical momentum that separates the regions of high and low energy is proportional to H, the square of the 
corresponding length being 2 (3A-i) TP- ^ e are S om g to see later that in the dc Sitter case this corresponds exactly 
to the Hubble horizon, the scale where the cosmological perturbations freeze out. Thus, the scale that separates the 
UV and IR limits in the " healthy extension" of HL gravity it is the same as the event horizon of de Sitter inflation. 
Beyond this length scale the evolution is similar to GR. 

We can write the second order Lagrangian for the perturbations using the constraints (|18p in the following way: 



S2S, 



dt- 



d 3 k 



■ a ( c 



fwipktyk + fw&kfyk - m lHl ~ rn^ipl - m^Vfe^fe) , 



(22) 



where the coefficients are given in Appendix 1. We can see that this action has two dynamical degrees of freedom 
with mixed kinetic terms. They cannot be combined into one degree of freedom, the Mukhanov-Sasaki variable, 
— £ = ip + -j-Scf), as is done in the case of cosmological perturbations in GR 0, H| or even in the "non-projectable" 



version of HL gravity with detailed balance [21 1 . 

We cannot diagonalize the entire action. Similar to what it is done in multi-field inflation, we are going to make 
a rotation in field space to find the true propagating and canonical variables of the problem (46|. For this we will 
diagonalize the kinetic term of the action. We can write the action in vector and matrix form as: 



S2S = x 



dt- 



d 3 k 



(2tt) 3 
8<Pk ipk 



Hk ipk 



C0 

2 C 4>4' 



2 C 0V' 



5fa 



( u 


f<M> \ 


(t) 


V f<H> 


U J 





4>k 



(23) 
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Diagonalizing the kinetic term which we call C, we find that the eigenvalues are given by 47| : 



d (k) Ai 



(3A-1) 



H 2 + 4> 2 



h (k) P 



(A-l) 



2(3A-1) 



± 



(3A-1T 

X 4 



£T 2 



(A-l) 
2(3A-i; 



-2r 



/i 2 (*) fc 4 



(A-l) 

2(3A- 1; 



n 2 



-2r 



(24) 



Given these eingenvalues we have to determine the rotation matrix S that diagonalizes C and then apply it to the 
entire action. Imposing the normalization condition, the eigenvectors of C can be determined and the rotation angle 
can be inferred from: 



tan(20) 



(25) 



We need to apply the diagonalization matrix in the entire action to write it in terms of the new variables. As a result 
of this algebra we find that, up to total derivative terms, we have that the action in terms of the new basis with 
diagonalized kinetic term is given by: 



S 2 S 



dt 



d 3 k 



ci</>3 + c 2 (f>l + F^fofo + F^fafa + Ffafafofa + Ffofafofa - M| 3 03 - M^J>\ - M^ 34>i (j> 3 (f> 



f2 ±2 



(26) 



where the new rotated variables are given by 



S*c 2 



(27) 



and the new coefficients are given in Appendix 2. 



III. COSMOLOGICAL PERTURBATIONS AND THE POWER SPECTRUM 



Because of the power-counting renormalizability, HL gravity can contain terms up to sixth order in spatial deriva- 
tives. The terms of highest order in momentum dominate the evolution in the UV. As pointed out in [35], HH this will 
lead to a scale invariant initial power spectrum of vacuum fluctuations. The reason is that the form of the state which 
minimizes the HL Hamiltonian is different from the one which minimizes the usual Hamiltonian of a Lorentz-invariant 
theory. The difference in Hamiltonians also affects the evolution of the state. We have to investigate the evolution 
of the perturbations that are generated inside the trans-Planck region until the IR region where we recover the GR 
evolution and verify how the new physics in the UV modifies the power spectrum. 

Another effect in the " healthy extension" is the presence of an extra degree of freedom that will only be dynamical 
in the UV region of the theory, inside the Hubble radius. The coupling of this extra mode to the adiabatic mode can 
also be a source of modifications of the power spectrum of the adiabatic fluctuation mode. In the following we will 
investigate these effects. 



A. Equations of Motion 



It is convenient to perform a change of variables to absorb the multiplicative factors in front the kinetic term and 
absorb the scale factor. The rescaled variables are: 

Uk = a^/Xi<f) 3 , v k = ay/X^4>4 . (28) 

Up to terms which are total derivatives, the action becomes: 



foS = X 2 / drj 



d 3 k 
'(2^)3 



-^0304° ' i 2 2 i 2 2 i f 2 



VA1A2 



(29) 
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where: 



2 3HF,p 3 



to,. 



-M£ v = —;= 
VA1A2 




With this, we can derive the equations of motion for both dynamical variables of the problem: 



■ ( 



Fa, 



i k -uj" u u k = -^{v k a 



VA1A2 



-A 



1 



F 



VAIA2 



+ Vk 



Mi 



VA1A2 



" 03 04 O 
VA1A2 



}■ 



(30) 

(31) 
(32) 

(33) 
(34) 



The evolution of these variables is of the form of coupled oscillators with time dependent mass, where the right side 
of the equations represents the interaction between the oscillators. 

After these calculations, some consistency checks must be done to see if the system describes the expected physics: 

• First consistency check: We need to check if the variable 03 is the Mukhanov-Sasaki variable, the one that we 
want to be the right canonical variable described by the perturbations in the IR. In order to check this, lets 
evaluate 03 e 04 : 

I,. n 1 

2 ( 



4>l = - (1 + cos 26)6(/>l + - (l-cos20)V>fc+ sin 29 Sfaipk, 
4>l = -(1 - cos20)<5<^ + -(1 + cos20)V>fe - sin20 6<t> k ip k , 



(35) 
(36) 



and using we can see that in the IR and when H 2 ^> ifio/x 2 we have: 

03 ~ 5<t> k + -J-1pk , 04 ~ — (50fe - tpk ■ 
00 00 



(37) 



We can see that the new variable 03 has the form of the usual Mukhanov-Sasaki variable in the IR and it will 
hence describe the right observable. 

• Second consistency check: Now, we need to verify that the extra variable 04 decouples in the IR, which means 
that its mass diverges in this limit, which then leads to the presence of a single propagating degree of freedom. 
Considering H 2 ^ 0q, we obtain that in the IR the eigenvalues are: 



and in the UV: 



Ai — > 



Ai 



1 

2? 



A- 



k 

a~H 



2 X 2 



A, 



,3A- 1 
' A- 1 



(38) 



(39) 



From (f3"3"l) , ([54]) and (|3"2"|) we can see that the mass terms of the equations are divided by the respective eigenvalue. 
So, from (|38|). we infer that A2 goes to zero in the IR when k — > and its mass goes to infinity. The extra 
degree of freedom 4 decouples in this limit, as required for the theory to have an IR behaviour compatible with 
observations. 



Having verified the consistency relations we need to evaluate the coefficients appearing in the equations of motion. 
We want to study the dynamics in the UV limit, as we are interested in investigating the trans-Planckian problem 
and the power spectrum at Hubble horizon exit. In this way, we will determine the coefficients in this regime, using 
the coefficients of Appendix 2 and (|38I39[) . 
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B. UV Limit 



In the UV, the most important contributions come from the highest orders in fc. Therefore we will only consider 
the highest order contributions in each coefficient: 

cx 0(l/k 2 ) ; F Mi cx 0(const.) ; M 2 3 cx 0(k 6 ) ; M| 304 cx 0(k 4 ) ; (40) 

F 04 cx 0(l/fc 2 ) ; F Mi oc 0(const.) ; M| 4 cx 0(fc 6 ) . (41) 

The contributions of F^, 3 and F^ are negligible in comparison with M 2 ^ and M 2 ^ , which are also present in ui 2 and 
in w 2 , respectively. Hence we can simplify the equations of motion in the UV by considering w 2 ~ — 2M? 3 + H + H 2 
and lj 2 ~ -2A/| 4 +H +H 2 : 

a 2 fc 6 /3 2 fc 4 2 2 , 2 \ , / 7 2 fc 4 



( + + cik 2 -U - % j u k = -v k aZ 2 + v k ( + S 3 a 2 J , (42) 

rfjL6 2;„4 \ /T 2 fc 4 \ 

+ ^ + 5 ' fc2 - W - ^ = -% a ° 2 + M M + fi 3« 2 J ■ (43) 

where in the interaction terms we also considered the constant terms (which are the leading ones) of i^ 3 ^> 4 and i*0 3 4 
since they are accompanied by a time derivative term. Without a careful analysis of the influence of these terms we 
cannot neglect them. The remaining coefficients are given in Appendix 3 in terms of the dimensionless coefficients. 
We can see that the terms have the expected scaling dimension since we know that [uj 2 ] s — 2. 

The coefficient functions of the Uk and Vk terms on the left hand side of the above equations determine the dispersion 
relations of the mode functions in the UV. The dominant term is the fc 6 term. These dispersion relations encode the 
change in the physics in the UV due to the presence of the higher space derivative terms in the HL action. These 
dispersion relations are of the Coley- Jacobson [3(| type and have been used in Q to study the trans-Planckian problem 
for cosmological perturbations in the context of toy models. 

As is obvious from the above equations, in the "healthy extension" of HL gravity we have to consider two effects 
in the evolution of fluctuations in the UV regime which are different from what is obtained using GR: firstly, the UV 
modification of the dispersion relation introduced by the UV complete action leads to a different form of the initial 
state and to a different evolution of the modes on trans-Planckian scales, and secondly there is an extra degree of 
freedom for scalar metric fluctuations which is active in the UV. Through the mode coupling, fluctuations in this 
second mode induce growing perturbations of the curvature fluctuations. If we imagine imposing initial conditions for 
all fluctuation modes at some initial time £j, then short wavelengths are subject to the presence of the source term for 
a longer period than l ong wavelengths. Hence, one might expect that mode mixing will lead to important deviations 
from scale- invariance [48j. 

The effects of modified dispersion relations on the spectrum of fluctuations in inflationary cosmology was studied 
in detail in Q (see also [13 for a review). It was found that if the dispersion relation is of the type for which the 
adiabatic evolution of the mode function is maintained, then no deviations from scale invariance of the perturbation 
spectrum results. The dispersion relation which wc obtain is in this class. Hence, we expect that deviations from 
scale-invariance of the spectrum of curvature fluctuations in an inflationary universe in HL gravity can only come 
from the mode mixing effect. 

After these preliminary discussions, we turn to the solution of the above system of coupled equations. In the linear 
approximation which we are using each Fourier mode evolves independently. The behaviour of the mode functions 
is different in the three regions of space-time indicated in Figure [1] Region I is the region where the fc 6 term in the 
dispersion relation dominates, i.e. where A -C l c . The wavelength l c where the k e term ceases to dominate is given 
by the mass parameters in the higher space derivative terms in the HL action. Based on cosmic ray constraints [38[ 
we know that l c must be smaller than the Planck length l p i . In principle there could be a small region of wavelengths 
where the fc 4 term in the dispersion relation dominates. However, this would require tuning of parameters. We will 
assume that we can neglect the term in the dispersion relation proportional to fc 4 , and this is justified if the conditions 
c 2 > IqP 2 (f° r Uk) and S 2 > IqC 2 (for Vk) are satisfied. Wc will hence assume that for A ^> l c the linear term in the 
dispersion relation dominates. Thus there is a region (Region II with l c <C A <C Ih-, Ih being the Hubble radius) 
where the dispersion relation is linear, as in standard GR theory. However, since the UV/IR scale is given by the 
Hubble radius, Ih, in Region II the extra degree of freedom is still dynamical and the interaction between the degrees 
of freedom must be considered. Region III is the IR domain where wavelengths are larger than the Hubble radius. In 
Region III where A> Ih the extra degree of freedom disappears. The perturbations will freeze out and lu 2 UI = a /a. 
Hence, in this region the evolution of curvature fluctuations is as in GR. Wc wish to calculate the power spectrum at 
the time 772 (fc) where the perturbation mode fc becomes larger than Ih- 



Region 1 Region II Region III 



FIG. 1: Space time diagram representing the perturbations modes in the UV regions: Region I, where the fc 6 term dominates 
and Region II, where Wk on k; and the IR Region III. In red we see the Planck length showing that the Region I must be 
sub-Planck. 



To construct the general solution to this system, we need to match the different solutions from different regions at 
the transition times. Considering an expansion of the universe of power-law inflation type, a(rj) = Zo?7 (note 
that we recover de Sitter inflation when b = 3), the matching must be performed at the times: 



1/2 u / x \ V2 



\rff\w = (jrj ' £, WI (1 - 6)/2 =(£) T ' (44) 

which is the transition between Regions I and II for the two modes, with i" = (2Tr)(a/c s ) 1 ' 2 and V c = (2n)(S/c s ) 1 ^ 2 . 
Considering the definition of the coefficients (|32|). we can see that: 

n 2\V* 




77 - 3g 8 + 2^ l pl = Cl p i . (45) 
7 H J 

From the constraints imposed by Gamma- Rays [38| l c and l c cannot be larger than the Planck length. We do not want 
them to be much smaller, either, otherwise the effective field theory description of HL gravity becomes questionable. 
Hence, we assume that the two lengths are a bit smaller but close to the Planck length, and hence C and C are of 
order one. 

The transition between Regions II and III takes place at the times 



l%l = ^^Zfc' = (46) 



with l u H = 47rc s Z |7 7 |( 3 - b )/ 2 /v / fo T ^ r T and l v H = ATr£ s l \r ] \ < - 3 ~ b '>/ 2 /Vb^ 



1. Initial Conditions 



Before solving this system of equation we need to determine the initial conditions. Since in inflationary cosmology 
initial classical fluctuations redshift, it is usually assumed that the initial conditions are given by quantum vacuum 
perturbations, and that the system that will evolve following the classical equations afterwards. In order to determine 
this initial state we will make use of the Hamiltonian formalism and canonically quantize the system. Since the 
Hamiltonian in HL gravity is different from that in GR on the UV scales which are relevant to us, the vacuum state 
(defined as the state which minimizes the Hamiltonian Q) will differ. 

From (|29|) . it follows that the conjugate momenta for the fluctuation modes are: 
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Hence we can write the Hamiltonian of the system for each mode: 

74 = p Uk u k +p Vk v' k - 5£ k = (u^ - uj 2 u u' k ) + (v' k 2 - uj 2 v v' k ) + M 2 v u k v k 



(48) 



The canonical quantization can now be done. We promote the variables to operators that satisfy the canonical 
commutation relations: 



[u>k(v),u k ' (t))] = \pu k (v),P u k' (v)] = , 

[Mv),P uk '(v)] = s (3) (k-k'), 



(49) 
(50) 



that come from [4>(x, t),p,p(y, t)} = i5^\x — y) (all the others are equal to zero). We have the same relations for v k - 
As the Region I is sub-Planck, the initial conditions ti arc given in the deep UV. In this limit, the term that 
dominates is the one proportional to k 6 and the expansion of the universe can be neglected. The interaction term 
is also neglected initially and our system can be considered as a free system in a Minkowsi space-time. With this 



Pu k 



2u k and p Vk 



2v t 



The fields of the non-interacting system can be expanded in terms of creation and annihilation operators in the 
standard way: 



u(x) = Vh 



d 3 P 

(2tt) 3 / 2 ' 
d 3 p 



d\e~ lkx 



(51) 
(52) 



where ur u = fc 6 a 2 /a 4 and uj^ 1 = fc 6 <5 2 /a 4 . We can write the Fourier components of the fields in terms of the annihilation 
an creation operators: 



i'k 




Puk 



Pvk 



2 



a k - a_ k 



(53) 
(54) 



The creation and annihilation operators obey the commutation relations [a^a^,] = <5^ 3 ^(fc — k ) and [ofc,o!/] = 

6^ 3 \k — k ), with all the others equal to zero. These commutation relations are valid if the following normalization 
condition holds: 



2i 



(55) 



which is the Wronskian of the classical solutions. This normalization allows us to fix the amplitude of itk (rj) and 
Vk (ti) such that it is compatible with the Heisenbcrg uncertainty principle. 

We can now determine the initial conditions for our modes. Because of the normalization conditions (|55l) . this task 
is reduced to the determination of the effective frequency at the initial time. As mentioned before, we choose our 
initial state to be the one that minimizes the energy. This is a good definition for the initial vacuum state since it 
generalizes how the vacuum state is usually defined to the case of a modified dispersion relation. As these fluctuations 
are initially oscillating like they would in a Minkowski space-time, our initial conditions prescription can be viewed 
as an application of the Einstein Equivalence Principle to our problem. Making the ansatz 



«k = nc(f?K 



7k 



vu = f k (?7)e 



Ok 



(56) 



we see that in order for these variables in this form to obey the commutation relations (f5"0"]) and the Wronskian 
condition (f5"5"j). we have that r k ^ k = 1 and f\^ k = 1. 
The energy of the system is: 



+ ^>k| 2 



.12 2 



.112-2 



(57) 
(58) 
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Considering the minimal quantum fluctuations allowed by the uncertainty principle, the minimum energy is E Utk (rji) 
ui^ and E v k {rji) = lu^ 1 which corresponds to: 



r k fe)=^- 1/2 , r*{ Vi )=u- 1 l\ (59) 
r'u(Vi)=0, ^(^ = 0. (60) 



With this, the initial conditions are: 



UkiVi) = , j. , = 1/27,3/2 , Vk(.Vi) = ; ffj ^ = rl/ 2t 3/2 ' ( 61 



a(r]i) v a(r?i) 



2. Solving the Equation of Motion 

Having determined the initial conditions we can now solve the equations of the system (|43[) in each region and 
then match them. The only degree of freedom that is propagating in the IR is u k and this is the observable one 
that we are interested in obtaining the power spectrum in order to compare with the results obtained in standard 
inflation. Because of this, we are only interested in calculating the solution of Uh. The perturbation v k acts as an 
" entropy" perturbation, analogous to what happens in multi- field inflation models [3l| , and will source the " adiabatic" 
perturbation u k in Regions I and II, decoupling in Region III. 

The equations of motion that we want to solve are: 

fa 2 k 6 

u k + — j- + c s fc2 ~— \ uk = S u , (63) 




» { 6 2 k 6 „o, 2 a"\ 
v k + — r + c s k - — )v k = S v . (64) 



T K . a -2 \ c _ „.' o , \ ^ K o „2 



where 

= -w fc aS 2 + v k (^-^ 2 — + s 3a^J , S v = -u k att 2 + u k \^2~ + ^3a 2 J , (65) 
are the interaction terms. 

The homogeneous equations of motion have the same form of the equations of motion studied in Q in the case of 
the Corky- Jacobson dispersion relation. In that work it was found that the evolution depends sensitively on the signs 
of the coefficients of the terms in the dispersion relation. In our case we know from the construction of the " healthy" 
HL action that the coefficients a, 5, c s and c s must be positive to avoid exponential, ghost and tachyonic instabilities. 
If we evaluate the " adiabaticity coefficient" for our dispersion relation in the UV region and in the case of de Sitter 
inflation, we obtain 



Q(k,v) 



^2 1,4 26 , „2/4 2 

(66) 



0,2/^26 _ ^4^2 



and it is hence clear that Q(k,rj) > 1 which implies that the adiabaticity condition is satisfied. Hence, we know 
that the state tracks the instantaneous vacuum state, and thus the homogeneous solutions approach those of GR 
for A > I pi. The homogeneous evolution alone does not lead to differences in the power spectrum of cosmological 
perturbations compared to what is obtained in standard inflation. 

On the other hand, in both Regions I and II these equations are going to have both a homogeneous and a particular 
solution, because these are the UV regions where the extra degree of freedom is dynamical and mode mixing occurs. 

Let us now turn to the effect of mode interactions in the two UV regions, and proceed to evaluate the analytic 
solutions of the equations of motion in the presence of the interaction in each region. For simplicity, we will choose 
a = 8 and c s = c s , in order to only have a common value of l c and In for both modes. 
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In Region I, the equation of motion reduces to: 



l Ik 



-uik = S u , 



a* 

v i k H —vi k = *>V . 



(67) 
(68) 



We will only show here the calculations for uik, since the solution for vik is identical, only changing the cocfhcicnts. 

We are going to solve this equation using the Green's function method with Born approximation. This means that 
our solution is going to be the general solution of the homogeneous equation Lu = 0, where L = (cP/dr] 2 ) + oj 2 , plus 
a particular solution of the inhomogeneous equation. The particular solution is constructed using the fundamental 
solution of LG(rj,a) = 8(rj — a), where G(j],a) is the Green's function. The solution is: 



Mv) = Ukiv) 



f u (a)G u (r),a)da 



(69) 



where u h is the homogeneous one that obeys the initial conditions (|62|) , and second term on the right hand side is the 
particular solution constructed from the Green's function. 
The homogeneous solution has two independent solutions: 



u h Ik {r 1 )=AM 1,2 Ji/2 b {z)+AM 1/2 J-i,2 b (z). 



(70) 



where z = g\rj\ b /b eg — ak 3 /l 2 . For the vjk equation we would have z — g\i]\ b /b where g = Sk 3 /! 2 .. To determine 
Ai and A 2 we use the initial conditions discussed earlier. Using the Wronskian of the Bcsscl functions, [J_„ J v -\ + 
Ju+\Jv\{z) = 2 sin[7r/(26)]/7rz, we obtain 



fJTV 



26sin(7r/(26)) 
gir 

' 26sin(7r/(2&)) 



IVi} 1 b UlkiVi) J-l/2b(Zi) 

g u' Ik (rn) Ji-±{zi) 

|??i| 1 ~ fc U Ik (Vi) Jl/2b(Zj) 

g u' Ik (rn) J i i(zi) 



(71) 
(72) 



unless 1/26 an integer, which we assume is not the case. 

In this regime, Zi is large and proportional to In/MVi) ^ 1- Thus, we can use the asymptotic expansion of the 
Bessel functions for large arguments, J u (x) — (2/ttx) 1 / 2 cos(x — vit /2 — 7r/4), and rewrite the coefficients as: 



■ f 7T 9\ 



1/2 



\Vi 



i(6-l)/2 



sin(7r/26)' 



±i 



\2b) m sin(7r/26) 



(73) 



where Xi(rj) = Zi(rj) + (7r/46) - (n/i) and yi(rj) = Zi{rj) - (7r/46) - (7r/4). 

We can rewrite the particular solution (in the case where we have two independent solutions for the homogeneous 
equation (f6"8")) ) as: 



Su(rj) = ui. k {rj) 



dr) W u 1 u 2 , k ('n )fu(r) ) - u 2 ,k(v) 



dr] W u 1 u Xt k{r) )Mv), 



(74) 



where W is the Wronskian of the solutions with index 1 and 2 which in Region I is W = —2AiA2bsm(n/2b)/n. The 
integration limits in Region I are from rji until the cutoff rf((k), which is the transition between Regions I and II. To 
write the interaction, S u we use the Born approximation where only the homogeneous solution of Vk is considered. 
The solution for v j k is the same as the one for u j k given in (|7D|) . 

We can now evaluate the particular solution. This is done in Appendix 4, where we considered the general case 
and the case of de Sitter inflation. Here we only discuss the second case, because in this limit we can evaluate the 
integrals exactly. 

As we can see from (|106p and (|107[) that the integrals only depend on the combination krj and that the integration 
limit kr]i = const. This shows the integrands for different modes are related by time translation, as is the late time 
cutoff. The cutoff corresponding to the transition between Regions I and II also obeys this symmetry. The only source 
of different evolution between different modes is the initial time which is fixed for all modes [49j . If the integrals are 
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dominated by earliest times, we would expect a the mode mixing to lead to a large deviation from scale-invariance, 
whereas if the integrals are dominated by the contribution close to the transition point between Regions I and II 
then scale-invariance will be preserved. We now have to check if the dependency on k which arises from the initial 
condition is strong enough to affect the power spectrum in an important way. 

The terms l/(krji) n in (|112[) . where n = 3, 2 and 6, respectively, are always very small since they depend on 
{X(r)i)/l c ) which is small since the initial length scale is much smaller than l c . Therefore these terms can be neglected 
and the solution is given by: 

Sui {m) = T jk^ {~B lce ~ 2lzl ( cos ( zi - i) cr + cos ( 2i - I) c + 

(cos (zi - I) c+ + cos (zt - ~) cj) I , (75) 

where cf = e ±M > ± (-l) 5 / 6 e ± ^ +l57r / 3 and cf = e ±«<+*W3(_i)7/6 ± e ±i Vi are oscillating terms in k. We can read off 
the important result that the solution maintains the overall scale invariance of the power spectrum since the overall 
amplitude is proportional to fc -1 / 2 . The phase z\ is also independent of scale. However, we have an oscillatory 
dependence on e ±JZi (since Zi = Zi(k)) coming from the intial integration limit (the initial time). Note that the 
oscillations have a frequency which is proportional to k 3 . 

Note that the overall amplitude of the inhomogeneous term is smaller than that of the homogeneous one, as one 
can check by inserting the expressions for the various constants. 



+ 



3tt7 2 1 



-2iz\ 



3V2 
16tt 2 



-2iz\ 



We now move on to a study of the solutions in Region 2 in which the equation of motion has a linear dispersion 
relation: 



l IIk 



c 2 s k 2 u IIk = S u , 



with a homogeneous solution given as a linear combination of two basis solutions 



(76) 



(77) 



As our homogeneous solution we take the one which matches with the full solution at the end of Region I (the effect 
of the mode mixing in Region II will yield the particular solution in Region II) . We need to match the solutions ui & 
and unk when the transition between the regions occur at 771, to determine the coefficients B\ and B 2 . This gives us: 



Bi=T 



uik(fli)e 



— ic s krji 



2sin(7r/26) 



Jh ,(l-b)/2 e ±ixi 



Bo 



ic s krji 



2sin(7r/26) r\i 



e T(ix/2b) e i Xl 



?—ic s krji 

2 

icskr)i 



Sui (771) 
c> 11/(771) 



—Sujir)! 
c s tz 



(78) 
(79) 



To determine the particular solution in this region, we use the same procedure as in Region I. The particular 
solution is calculated in Appendix 4. Again, considering the case of dc Sitter inflation, it follows from (|114[) that 
both of the integration limits are constant, krji = (cs/oO 1 ^ 2 ^ — const, and krj2 = V%/c s = const. .The integration 
range will hence be independent of k. On the other hand, a dependence on k enters via the homogeneous solution 
in the source term. The k-dcpendcncc of the particular solution comes completely from the coefficients B\ and B 2 
(the coefficients analogous to the above B\ and B 2 which appear in the homogeneous solution for Vk) which are both 
proportional to /c -1 / 2 , indicating that the overall amplitude of the inhomogeneous term will respect scale invariance. 
However, like in the contribution to the inhomogeneous term in Region I, we still have oscillation. 

Evaluating the 9 integrals present in the solution (|114p , and using the asymptotic expansion of the exponential 
integral function for large arguments, Ei(z) ~ e _z / ' z + ... [39), pg. 231] we consider the dominant contribution from 
each term. Considering that the contributions of the particular solution in Region I and its derivative in the coefficients 
B\ and B 2 , given in (|79[) . are smaller that of the homogeneous one, we can neglect them. With this approximation, 
the particular solution becomes: 



5u n (n 2 ) = ±- 



J_Ih 

2tt 2 L 



7r 2 7 2 V2 2 / 



i(2V2) 



(80) 
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where 



a = sin(v / 2(^Ac)+yi)-sin(v / 2(^//c) + a;i±7r/6), 
a = cos(V2(l H /l c ) + yi)-cos{V2(l H /lc)+x 1 ±Tr/6) 
A = sin yz (1 + Aijjii c ) + m -sin 



(81) 



\/2 (1 + Hh/Ic) + yi\ - sin \V2 (1 + 3l H /l c ) +xi± tt/g] , 



are constants containing terms that depend on z\ which is independent of k. 

Comparing term by term, this particular solution dominates over the one from Region I. The particular solution 
is oscillatory in k (with a frequency which is proportional to k 3 ) via the overall coefficient e ±lXi . Comparing the 
amplitude of the particular solution with that of the homogeneous one, we see that the overall amplitudes (outside 
the square bracket of (|80|) ) are the same. The relative amplitudes of the oscillatory and constant terms in the mode 
functions is given by the coefficients inside the square bracket of (|8T?|) . Taking the new scale appearing in the HL 
Lagrangian to be comparable to the Planck scale we find that the amplitudes of the first two terms are of the order 
one, not suppressed by factors of H/m p i as the oscillations obtained in (ll| are. The amplitude of the third term is 
suppressed by the inflationary slow-roll parameter. 



Since 



V2 (Jc 



3 \X(Vi)J \KVi) 



( In 



(82) 

we see that the frequency of oscillation is set by the initial time and is very large in units of ^fej (since the second 
factor on the right hand side of (|52"j) is much larger than unity). Hence, these oscillations will be very hard to detect 
observationally given a finite frequency resolution of an experiment. 

To evaluate the power spectrum of curvature fluctuations we need to have the solution in Region III and match it 
with the full solution in Region II at Ijj crossing. 

In Region III, we have the equation: 



u nik H u nik = 



(83) 



with solution Umk = C a(rj). Now, we need to match the solutions when the mode exits the Hubble radius at 772 to 
determine C . At 772 the solution in Region II is given by: 

uu k fa) = T ^777^ e ±iXi (e^^-^h + e-^'-^h*) 



1 



Suifa) -8ujfa) ) + Sun fa) ■ 

Co K 



(84) 



where h = e iyi — e ' l ( x i ±7T / 2 ) and h* it complex conjugate. Thus, the coefficient is C = uufa)/a. 



3. Power Spectrum 



Finally we can evaluate the power spectrum of the fa, the Mukhanov-Sasaki variable . This is the only propagating 
mode in the IR and will be the " adiabatic" perturbation produced during inflation. The dimensionlcss power spectrum 
will be calculated at Hubble horizon crossing. Afterwards the perturbations freeze until re-enter the Hubble horizon 
later in the evolution of the universe. The power spectrum is given by: 



fc 3 P = k 3 \Cf 



k 3 r). 



I 2 

47T 2 fc 



unfa)Y 



J/2 



k l/2 



Jc s (V2-ri2)f l _|_ g- ic s (V2-V2)fti 



3 ic s (j)2—m) 



-ic s (r?2— J72) 



Suifa 



-Sujfa 



Sunfa) 



(85) 



We can see from this solution that the amplitude of the power spectrum will be scale invariant since all terms are 
proportional to fc 1 / 2 , as we can see in (|75|) and (|80|) . As already pointed out, we will have an oscillation in this 
power spectrum coming from terms that depend on Zi and its mixing with the other terms in the sum. The relative 
amplitude of the oscillations is of order one, show a characteristic A:-dcpcndcnt frequency, but are too rapid to be 
detectable by an experiment with finite frequency resolution. 
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IV. CONCLUSION 



We have studied the generation and early evolution of cosmological fluctuations in the "healthy extension" of 
Hofava-Lifshitz gravity, assuming that there is matter which leads to an inflationary background cosmology. Since 
HL gravity can be viewed as a consistent UV completion of gravity, it being a power-counting rcnormalizablc model 
for quantum gravity, inflation in the context of HL gravity provides an interesting background to study the trans- 
Planckian problem for inflationary cosmological perturbations. 

There are two reasons which could lead to large deviations from the usual predictions of standard inflationary 
cosmology: firstly, the dispersion relation on trans-Planckian scales is non-standard, and secondly there is a second 
scalar metric fluctuation mode which acts as a source term in the evolution of the curvature fluctuations. The time 
interval during which mode mixing occurs depends on the wavenumber k of the fluctuation mode. 

The result of our analysis is that the overall scale-invariance of cosmological perturbations is maintained in HL 
inflation. On the other hand, there are oscillations in the amplitude of the spectrum as a function of k whose relative 
amplitude is of order one, and whose frequency shows a characteristic scaling proportional to A; 3 , the power being 
determined by the power of the highest power of the spatial derivative term in the HL action. The frequency of 
oscillation is so high, however, that it is unlikely that an experiment with finite frequency resolution will be able to 
detect these oscillations. 

The reason why scale-invariance of the spectrum of cosmological perturbations is maintained in HL inflation is that 
the dispersion relation for the fluctuation modes satisfies the adiabaticity condition, and that the non-trivial initial 
vacuum state evolves into the usual one once the wavelength drops to sub-Planckian values. The oscillations in the 
spectrum come from the particular solution of the mode evolution equation (induced by mode mixing). Physically, 
the oscillations are a consequence of the fact that setting initial conditions for all Fourier modes of the fluctuation 
variable at a fixed initial time breaks the time-translation invariance of the mode evolution. 
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Appendix 1: Coefficients of the Second Order Action 

The coefficients of the second order Lagrangian of the Healthy extension of Horava-Lifshitz gravity are 



d(k) ClP = 2(3A-l)^ r + (A-l)/ 1 (fc)^ J 



d (k) fa 



d (k) c v , = 2 (3A - 1) 



% + 2/ 1 (fc)P 



d{k)c^ = 4:{3X-1)H^ 



d (k) U = - (3A - 1) - (A - 1) [V ,4, + Vi (0o) k 



X 



41 90 

X 4 



-24 (3A - 1) HA + 12 (3A - l) 2 H 3 - 6A (3A - 1) - 12 (3A - 1) ?± U ) H . 

X X 



d (k) f H , = 6 (A - 1) %A - 3 (3A - 1) (3A + 1) ^H 2 - § (A - 1) 4 + 



3(A-1 
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X 



X 



X 



ff-2(A-l)[3/ 1 (fc)+2/ 2 (fc)]^P, 



d{h)f H = 4 (3A - 1) ^ OM H - (3A - 1) || - 2 (3A - 1) f x (fc) M + 4 (3A - 1) Va (0 j ^ 



(86) 

(87) 
(88) 
(89) 
(90) 
(91) 



Appendix 2: Rotated Coefficients 



The coefficients of the rotated second order Lagrangian are 



i(l + cos20)/ + 1(1 - cos20)/ v , + S -^(f 4 



1(1 -cos 20)/^ + \i l + cos29)U 



2 

sin 26 



sin 2$ 1 1 
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1 , , , 1 . , 9 sin 20 9 
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where: 



sin 20 



cos 20 



.2 

-01/, 



^ (C^ - C ) 2 + ^ (c^ - C0) 2 + C 2 

The general form of the c x coefficients it quite involved and will not be written down here. 



Appendix 3: Coefficients of the UV EOM 

The coefficients of the UV equation of motion are given by: 



2 A v 6 

a =4-r, 
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Appendix 4 : Particular Solutions 



Region I: 
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We are going to evaluate the particular solution, Suj =Ti + T 2 . Making the change of variable rj = krj, we have 
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case of de Sitter inflation, this particular solution reduces 
i /7T\ 3 / 2 77 1 / 2 . . 



where z = ct/c 3 b rj b /bl^. In the 
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+S3Z 2 
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■%(fc) -/)7»(fe) 
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(107) 



where z = a?7 3 /3Zo- 

We can now evaluate the 12 integrals that we have in this solution. We are only computing them for de Sitter 
inflation. The solution of these integrals are polynomials multiplied by generalized hypergeometric functions. This 
functions can be written in terms of Meijer functions, a general function that contains most of the special functions 
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as a special case and that it is defined as a path integral in the complex plane [40l . |4l| . For example, the integral: 

f*7i 

-2 



ii = ts = 
fj 3 



drjT] J_i/26 (z) Ji/2fc (z) 

iF 3 (1/2, 1/2; 5/6, 7/6, 9/6; -z 2 



Wl,2 / ,2 

6^ 2 < 4 



1/2,1/2, 
0,1/6,-1/6,1/2 



(108) 



(109) 



All integrals give the generalized hypergeometric function with p = 2 and q = p + 1 = 3 which is a special case of the 
Meijer function with m = 1, n = 2 = p and q = 4, except for the last integral in (|107[) (and its counterpart in (|106p ) 
that gives a hypergeometric function with p — 3 and q = 4 that it is a special case of a Meijer function with m = 1, 
n = 3 = p and g = 5. 

To be able to work with these results we make use of the asymptotic expansion of the Meijer function for large z. 
For the Meijer functions present here, we use Theorem 1.8.5, the case (1.8.13) in [4l|, verifying that the functions 
obtained obey all the conditions imposed in the theorem: 



G 



n,m / 2 



■}••••> ^m? 3 * • • ) 

&1, . . . , &„, &71+1, . . . , b q 



D™(\)H M [z 2 exp^(^-2A)]] , 



(110) 



where, jj* = q — m — n, v* = — p + m + n and arg(z 2 ) = 0. In all of our cases fi* = 1 = /i* and Dp l q l (\) = 



A rn,n = ^y* exp ^ fl . _ ^] =m+1 b 3 

be written as a series: 



as follows from (1.7.23) and (1-7.2) of [41| . The function ilp, g can 



ffp, g [a;] = exp (p - g)x"-« 



(2^ 



(g-p)V2 



^ +0(.T 



(111) 



with p* = [EUn+i h - E^i a ] + (p-q+ l)/2j /(g - p). 

This asymptotic expansion was applied in all of the results and, after some algebra, the particular solution simplifies 



to: 
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COS I 01 I cj + COS I Z\ 
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TV 



(112) 



Region II: 

To determine the particular solution, Sun = Ti + T 2 , we make the change of variable fj = kr/, and get: 
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(113) 



where B\ and B 2 are the coefficients of the homogeneous equation in Region II for Vk- They are equal to (|79|) . since 
the coefficients of the equations of motions of u and v were set to be the same, for simplicity. 
For de Sitter inflation, the solution is: 
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